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Abstract 

The correlation function of two dimensional Ising model with the nearest neigh- 
bours interaction on the finite size lattice with the periodical boundary conditions is 
derived. The expressions similar to the form factor representation are obtained both 
for the ferromagnetic and paramagnetic regions of coupling constant. The peculiarities 
caused by the finite size of the system are discussed. The asymptotic behaviour of the 
corresponding quantities as functions of size is calculated. Some generalization of the 
obtained result is conjectured. 
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1 Introduction 



Since the outstanding result of MontroU, Potts, Ward [1] there appeared many papers de- 
voted to the problem of the spin-spin correlation function in the two dimensional Ising model 
(IM) with the nearest- neighbour interaction on the infinite size lattice (see e.g. refs. in [2]). 
The basic IM theory results are collected in the significant book [3]. I can not, of course, to 
give here any comprehensive list of IM references. Only very few of them which were straight 
instructive for the preparation of the present paper will be referred to. I think, in particular, 
that the results and ideas of the paper [4] were very fruitful for the further development of 
the IM theory. The study of the scahng hmit of the IM (see e.g. [4]-[7]) is of interest from 
the quantum field theory point of view. 

The so-called form factor representation [7] for the two dimensional IM correlation func- 
tion appears to be a crown of this activity: 

((7(ri)(7(r2)) = const ^c/„(r), r = |ri-r2|, (1.1) 

n 
oo 

Fn[q] = fl(^^): a;, = uj{q,) = ^m^ + ql (1.3) 

In the ferromagnetic region of the coupling parameter the summation in (1.1) is extended 
over even n, in the paramagnetic — over odd. It is worth of noting that this representation 
was first evaluated [7] in the framework of the iS-matrix approach [8] and then [6] by means 
of straightforward IM solution. The discovery of the form factor representation was very 
fruitful for the progress of exactly integrable quantum field theories [9]. 

The advantage of representation (1.1) consists in that the dynamical aspects of the system 
are separated from the kinematical ones: form factors squared are integrated over the phase 
volumes of the n-particle intermediate states. The form factor representation (1.1) clarifies 
the dynamics of the model but does not answer on the questions about the spectrum, sort 
and statistics of the particles that form the intermediate and asymptotic states of the system 
under consideration. The analysis of the model at finite (nonzero) temperature or in the 
finite volume is necessary for this purpose. One can observe an activity in this field in last 
years [10]-[16]. These works show in particular that the problem is complicated enough: the 
authors of refs. [15], [16], for example, call the conjectures of [12], [14] in question. I think 
that a simple exactly solvable lattice model example would be very useful for the business. 
Formally the finite temperature (in quantum field theory sense) means the finite size along 
the temporal axis and periodical boundary condition for boson fields or antiperiodical — for 
fermion fields. Meanwhile there is no representation which is analogous to the form factor 
one for the correlation function on the finite size lattice even for the Ising model. 

This work, I hope, makes up a deficiency. I have calculated the IM spin-spin correlator 
on the lattice wrapped on a cylinder by the use of the classic methods of IM theory [3] 
adapted properly to the case of the finite sizes. The solution is expressed in the form similar 
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to the form factor expansion (1.1)-(1.3). For reader's convenience I write the result just in 
the Introduction. If one considers it obvious he would be free of cumbersome mathematical 
transformations. So, the expression obtained in this work is 



(1.4) 



n 



9n{r) 




(1.5) 



n\{N) 



n 



Fn[q\ = n 



n 



r sin((gi-gi)/2) 
sinh((7, + 7,0/2) ■ 



(1.6) 



More precisely, the representation (1.4)-(1.6) corresponds to the correlation of spins posed 
on the line parallel to the axis of a cylinder with N sites on its circumference. The values 
^, ^T, A, 7j, rii, in eqs. (1.4)-(1.6) are defined further in Sections 2, 3, 4. Three of them — 

A and rji — are specific cylinder circumference dependent values: In ^7-, A~^, r]i vanish 
if tends to infinity. The appearance of the summation in (1.5) instead of the integration 
over phase volume is natural consequence of the size finiteness. It is widely believed that the 
underlying IM field is the fermion one: for instance the IM partition function is exactly the 
same as in the free Majorana fermion system. So, the boson spectrum of quasimomentum in 
Brillouin zone (this is denoted by the upper index in sum (1.5)) is surprising . Nevertheless 
it follows unambiguously from the calculations. It takes attention also that the lattice form 
factor (1.6) does not depends explicitly on the cylinder circumference. Note that similar 
expression for at even n was found by authors of [6] for the infinite lattice case. 

In Section 2 the model is formulated and the brief evaluation of the Toeplitz determi- 
nant representation for the IM correlation function is given. The corresponding Toeplitz 
matrix allowing for the finite size lattice is calculated. In Section 3 the lattice form factor 
representation for the correlation function is derived for the ferromagnetic domain of the 
coupling parameter and for the paramagnetic domain — in Section 4. In Conclusion I dis- 
cuss the relationship between the IM correlation function on a cylinder and quantum field 
Green function at finite temperature or in finite volume. I also comment on some prospect 
of generalization of the obtained results. In Appendix the factorised representations for 
the Toeplitz determinants of special form which are used for obtaining lattice form factor 
expansions are evaluated. 



The Ising model with the nearest neighbour interaction on the square M x N lattice (see. 
Fig. 1) is defined by the hamiltonian H[a] 



2 The Model 



PH[a] = Kj2a{r){V, + V,)a(r), K = (3 J, 



T 
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where J is the couphng parameter, (3 is the inverse temperature; two dimensional vector 
r = {x, y) numerates the lattice sites: x = 1, 2, . . . , M, y — 1, 2, . . . , N] Vi, Vj, are the 
one step shift operators 



Y 

N 



V^(7(x, y) = a{x + 1, y), Vya{x, y) = a{x, y + 1). 
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Fig. 1. The positions of the correlating Ising spins on the lattice. 



The partition function of the model is 

Z = ^e-^^M, (2.1) 
M 

the pair correlation function is 

{a{v,)a{v2)) = Z-^^e-^^M<j(rO<7(r2). (2.2) 

M 

For the lattice with the periodical boundary conditions (wrapped on a torus) the par- 
tition function (2.1) can be expressed [17], [18] through the combination of the four Gauss 
Grassmann functional integrals with different boundary conditions 

Z = (2cosh2K)^^i(Q(^/) + + - Q^"")), (2.3) 

where 
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The action S[ip] is the anticommuting quadratic form 

r,r' i',p=l 



D 



V 



l + tVa: 


1 


1 


\ 


1 - tv_^ 


-1 


1 




-1 1 









-1 -1 







/ 



t — tanh K. 



(2.4) 



(2.5) 



The indices (/, b) in eq. (2.3) shows the type of the boundary conditions for the shift 
operators in the matrix (2.5). 

7" - antiperiodical: (V^^))^ = -1, (Vj^))^ = -1; 

the corresponding quasimomentum runs over halfinteger values in Brillouin zone 
(— TT, tt) (in 27r/M unites for and 27r/N - for py) 

"6" - periodical: (V^**))^ = 1, (Vf )^ = 1; 

the corresponding components of quasimomentum runs over integer values. 

The correlation function (2.2) can be also presented in terms of the Grassmann functional 
integrals with the Gauss distribution [18] 



where Qd is the functional integral 



Q(ff) + Q{bf) + Qifb) _ Q(bb) ' 



(2.6) 



(2.7) 



The action with the defect denoted by Sd[ip\ in (2.7) differs from (2.4) in that the parameter 
t is replaced by t'^ on the links along the path which connects the sites ri and Y2. 

I consider the special case when the sites ri and V2 are situated on the line parallel to 
the horizontal axis as is shown on Fig. 1: 

Ti = {x,y), r2^{x + r,y), |ri - = r. 

The action with the defect in this case has the form 

^ S[^] + AS[^|;], 

where 

r 

AS[ip] = (r^ -t)J2^\x + x'- 1, y)ip^{x + x', y). 

x'=l 

The ratio in the r.h.s. of (2.6) simplifies 



Qiff) 

(<7(rO<7(r2)) = rg^[l + o(e-^/^)], 
5 



(2.8) 



if at least one of the lattice sizes is much larger than other, i.e. torus degenerates to cylinder: 
in our case M/N ^ 1. Note, that only the terms with the antiperiodical boundary conditions 
survive. 

After some manipulations with the functional integrals the ratio in the r.h.s. of (2.8) can 
be transformed to the functional integral over Grassmann field which is defined just on the 
line connecting the sites ri and r2 

{a{v,)a{v,)) = J d\x]e^^'^^'^ = det A, (2.9) 

r-l 

d[x] = n (^xl dxl, ix^Ax^) = J2 xlAx,x'xl', 



x=0 



^ ^ J_^if) e^^-("-"'^[2t(l + e) t^)(e-P- + t^e-^P-)] 

MN^ {1 + tY -'2t{l-t^) {cos + cos Py) ' ^ ' ' 

x,x' — 0, 1, . . . , r — 1. 

I recall that the index (/) points the antiperiodicity: the momentum's components run over 
halfinteger values - fermion spectrum 

E^'^Px) = E-(|(^ + E'^'viPy) = E.(|(/ + i) 

1=1 ^ Py 1=1 ^-^^ 

Now let us rewrite the expression (2.10) in appropriate for the Toeplitz determinant 
theory manner. Prom the tabulated formulas [20] 

n^^^2(cosh7-cosp^) = e^^(l + e-^^)2, I]^'^2(cosh7 - cosp^) = e^^(l - e-^^)^, (2.11) 

Py Py 

it follows in particular 

1 yif) 1 _ tanh(iV7/2) 1 y{b) 1 _ coth(iV7/2) 

N'^ cosh 7 — cos sinh7 ' cosh 7 — cos sinh7 

At the cylinder case M ^ 00, N — const the summation over p^ turns into the integral 

\z\=l 

Summing up (2.10) over py by use of (2.12) and accounting for (2.13) one can obtain for the 
matrix A^^^' 

= hf f (2-") 

where the kernel A[z) is 



A{z) 



a-e-7^z)(l-e-7°z-i) _ 
^(l-e-^o^)(l-e-T^O ' ^^-^^^ 
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and 

r(z)=tanh[(7V^(p))/2], (2.16) 

7o=|2X + lni|, 7^ = 2X-lnf. (2.17) 
The function 7(p) in (2.16) is defined by the following equation 

4sinh2^ = ^2^4sin2^, (2.18) 

where 

ii^ = 2(Vsinh2is: - 1/Vsinh 2^)^. (2.19) 
It follows from eqs. (2.18), (2.19) in particular 

7(0) = 7o, 7(77) =7^, cosh77r - cosh7o = 2. 

The function T{z) (2.16) tends to unity at the thermodynamic limit N ^ 00 (more 
precise N/j, :$> 1), and the kernel A{z) (2.15) turns into the classic IM theory expression. In 
this limit the value K = is the critical point: the specific heat and correlation length 
diverge. At the finite the corresponding singularities are smoothed and there is no phase 
transition. Nevertheless also in this case the value K = Kc is notable. The calculation 
technique based on Wiener-Hopf sum equation "distinguishes" the ferro- and paramagnetic 
regions and the expression for the correlation function at K > Kc differs from that at 
K > Kc- Just in this sense I shall mean further the words "critical point" , "phase" and so 
on. Note that the following equalities take place at the critical point 

sinh 2Kc = 1, A* = 0, 70 = 0, 7^ = 2 ln(v^ + 1). 

3 The Ferromagnetic region: K > Kc 

The dependence of the correlation function (2.9) 

((7(ri)a(r2))=detA('-) (3.1) 

on the distance r — |ri — enters through the matrix dimension denoted by upper index 
(r) in the r.h.s. of (3.1). The dependence on the cylinder circumference N comes into the 
function T{z) (2.16) and the dependence on the coupling parameter K is defined by the 
parameter n (2.19) in the eq. (2.18) for the function 7(p). 

The kernel (2.14) of the matrix (2.9) A^^^ has to be represented in the factorized form 
according to the method of the determinant calculation which we use (see (A6), (A7) in the 
Appendix) , 

A{z)^P{z)Q{z-'), 

where the functions P{z) and Q{z) are analytic ones inside the \z\ = 1 circle. It can be 
formally done by use of the projection operators (A8) that are defined in the Appendix 

2\nP{z) = \n(l-e-'^-z) -\n{l-e-^°z) +2VlnT{z), 
2\nQ{z-^) = \n{l-e-'^°z-^) -\n{l-e-'^^z-^) + 2VlnT{z). 
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Otherwise accounting for the exphcit form (2.16) one can express T{z) through the ratio of 
the products with help of eqs. (2.11) 

n''^'' [cosh 7(j9) — cos q] 

T\z) = . (3.2) 

n [cosh7(p) — cosg] 

It follows from (2.18) 

cosh7(p) — cosg = 2 + — — cosg — cosp = — e~'^^^^2;)(l — e"'^'^'^^ z~^) , 

so eq. (3.2) turns into the form 

j]We7{9)(i _ e-T(«)^)(l - c-T('')^-i) 

T'^{z\ = - (3 3) 

^' n(/)e7(«)(l -e-7(«)^)(l -e-7(«)^-i)' ^'^ 

Q 

Owing to the momentum q runs the entire Brillouin zone all the factors of the r.h.s. (3.3) 
are double except for the values q = and g = tt in the nominator. So the functions A{z), 
P{z) and Q{z~^) have no branch points rather the simple poles and zeroes: 

PW=exp(^-j )- nW(l-e^^(»).-) - 

0<g<7r 0<q<7r 



According to (A25), (A28) 



detAW^e'^W^^/aKO- (3-5) 

1=0 



Let us write the function h{r) as the sum of three terms 

/i(r) = -r/A + ln^ + ln^T- (3.6) 
It follows from the eqs. (A23) and (2.15) 

A"^ = -[lnP(0) + lnQ(0)] = --^ i — In T{z) ^ - [ dp In coth{N-f{p)/ 2). (3.7) 

2m J z TT J 

\z\=l 

The parameter A has the sense of the coherence length. Its asymptotic behaviour at large 
N is 

A:^e^^°./-^^. (3.8) 
V2sinh7o ^ ^ 
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One can see from (3.8) that the coherence length grows very rapidly when the cylinder 
circumference increases and 70 7^ 0. The lattice can be considered as infinite at ^ 
max[ln(A^7o), ln(r7o)]. It follows also from (A23) that 



In(^^r) = — <fdz\nQ{z-^)-^\nP{z 



2m J 

\z\=l 



dz 



(27rz; 



d{ZiZ2 



\Z2\>\Z1\ 



\nA{zi) \nA{z2) 

{Z2 - ZiY 



(3.9) 



or separately for the In ^ and In 



ln£ = - In 

4 



In^T = 



{2m) 



sinh 7o ■ sinh 7^ 
sinh2((7o + 7^)/2) 

d{ZiZ2 



22|>|^l| 



Z2 - Zif 



\nT{zi)\nT{z2). 



(3.10) 
(3.11) 



The r.h.s. of the eq. (3.10) does not depend on N and is singular at the critical point 



In ^ ~ - In at // — > 0. 



(3.12) 



For the calculation In^^ let us return to the momentum variables and integrate (3.11) by 
parts. The result is 



or 



In^T = 



In^T = 



dp dq '-f' {p)'~f' (q) 



271"^ J sinh(iV7(p))sinh(A^7(g)) 



iV2 



d'Jp d'Jq 



In 



In 



sin((p + q)/2) 



27r2 7 sinh(A^7(p))sinh(iV7(g)) 
Differentiating eq. (3.13) with respect to /i one can obtain 

d . . 1 



sm{{p-q)/2) 
sin((p + q)/2) 



sm{{p-q)/2) 



where 



ln6 



N f 
^^^^2^ J si 



[Ci(/X) -C2(/X)], 



Cl 



d'-f 



I cosh 77r — cosh 7 



C2(/^) 



N 



sinh(A^7) y cosh 7 — cosh 70 ' 



(3.13) 



(3.14) 



(3.15) 



d'y 



' cosh 7 — cosh 7o 



27r J sinh(A^7) y cosh77r — cosh 7 
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Deriving eqs. (3.14), (3.15) one has to account for that the integrand of (3.13) vanishes on 
the boundaries of the integration region. Note also that the derivatives of p and q are 

dp 1 dq 1 

2sinp' dii"^ 
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2 sin g ' 



that follows from (2.18). For A^7o ^ 1 the asymptotics of ci(//) and 02(11) are 



'^^^^"cosh('iV7o)V^^ri^' ^^(^)"4cosh(iV7o)V^- ^^"'^^ 

Collecting (3.14)-(3.16) one can obtain 

00 

In^T = / di^i^[cl{i^) - clip)] - -e"'"^^". 

So, we see that — > 1 if iV — >• 00. 

The integral (3.13) diverges at the critical point. This means that In^r as a function of 
II has a singularity. It easy to obtain at // — > 

^In^T 1 



and consequently 

In^r = — - In/x + const. (3-17) 

Therefore, one can see from (3.12), (3.17) that the function In(^^T) has no singularity at 
/X = if iV is finite. 

Let now evaluate the expansion in the r.h.s. of (3.5) 

00 

GF{r) = Y^g2iir). (3.18) 

1=0 

The coefficients g2i{r) are expressed through the 2Z-multiple contour integrals (A30) as it is 
shown in the Appendix 

21 i-i I I 

r_lV . UidZiZDU n [iz2i-l- Z2j-l)\z2i- Z2j)^] UW{z2i-l) 

m\{2m)'^J^^ n n (1 - ^2^-1^2,)^ n w{z2^) 

i=l j=l i=l 

The integral (3.19) is defined by the concrete form of the function W{z) = P{z)/Q{z~^): in 
our case this is the eq. (3.4). The integrand in the r.h.s. of (3.19) is analytic at \zi\ < 1 
except for the singularities which the functions W{zi) and W~^{z^^) possess. To separate 
explicitly the singularities inside the circle j^j = 1 these functions ought to be rewritten in 
the following way 



Wiz)^ ^^^^ 



l_e-7.^)(l- 6-7.^-1) P2(^) ^^^^^ 



Q{z-^) \| (1 -e-Toz)(l -e-To^-1) T(z 

W-\z-') 



(l_e-7o^)(l_e-7o^-i)g2(^) ^^^^^^ 



P(z-i) \ (1 - e-T-^)(l - e-T-^-i) T(z) 
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The analytic at \z\ < 1 functions Pt{z) and Qt{z) are defined by tfie equations 

InPT(^) = VhiT{z) = — /^^lntanh(7V7(p')/2), (3.22) 

\z\=l 

]iiQt{z-^) = 7'lnr(^) = — : i -lntanh(Ar7(p')/2), (3.23) 

27ri J z — z 

\z\=l 

z' = e-'P'. 

Note that it follows from (3.7), (3.22), (3.23) 

InPT(O) = -1/A, IngT(O) = 0, Pt{z) = Pt(0) ■ Qt{z). 
With allowance of these eqs. one can obtain instead of (3.20), (3.21) 

W{z) = e(^°-^'')/^-^/^(cosh7. - cosp)Q'Azf-^^^^^PW^, (3-24) 

smh 7(p) 

W-\z-') = e(^--^»)/^(cosh7o - cosp)g^(^) '"^^/^^l^y^\ (3.25) 

sinh7(p) 

coth(7V7(p)/2) _ 1 (b) 1 

sinh7(p) cosh 7(g) — cos p 

It is seen from (3.24)-(3.26) that the singularities inside the integration contours of (3.19) 
are the simple poles at the points 

^W=e-^(^«), qik) = ^, k = l,...,N, 
corresponding to the boson spectrum of the quasimomentum q. Accounting for 

1 rdz v{z) _ ^(e-T^^)) 



2m J z cosh 7(gf) — cos p sinh7(g)' 

\z\=l 

where v[z) is analytic at \z\ < 1, the integral (3.19) is fulfilled by residua. The result is 

...W=(^E"n(-i^)&.b]. (3.27) 

where 

G2i[q] = Cf^n n [(e-^---e-^^-0'(e-^"-e-^^0']x 

i=l j=i+l 

I 

n [(1 + cosg2i-i)(l - cosq2i)] 
XT — ' (3-28) 

n n n (1 - e-'^-'iy 

i=l i=l j=l 
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7^ = 7(?0, 77i = r^(?0 = -lng|(e-^(^*)), (3.29) 

C*;^' = — is the binomial coefficient. 

' IU\ 

Using the equahties 

^{e^^' — e~'''-'')(l — e^'^^'^) = cosh7j — cosh7j = cosg^ — cosgj, 
the eq. (3.28) can be transformed to 

Q2i[q] - -ji , 

nsinh2((7, + 7,)/2) 

i<j 

where 

i-i I 

V2i[q] = 2-2''cf n n [(cos?2i-i -cos?2j-i)^(cosg2i -cosg2j)1 X 

i=l j=i+l 

I 

X n(l + cosg2i-i)(l -cosg2i). (3.30) 

i=l 

The expression (3.30) is not symmetric with respect to changes the summation variables Qi 
with even subscripts by that with odd subscripts. The summation extracts the symmetric 
part of the function V2i[q] which is denoted by l2f™[g] 



V, 



'-'I p 



where the sum over all permutations of q2i-i and q2j (theirs number is just Cf'') is denoted 
by J2 ■ The equality 

V2r[q] = Ur^'iq], (3.31) 

is the crucial for the final result. Here the even (with respect to Qi <-> — g^) part of the 
product 

U2i[q] = f[sin'^'^'-'^^ 

i<j 

is denoted by U^r^k] 



21 



[o-=±l] i<j 
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Consequently the function V2i[q] under summation (3.27) can be replaced by the function 
U2i[q]- This makes it possible to present the coefficient g2i{r) in the form similar to that of 
eq. (1.2) 

n = 21, 



-n/A ^^^^ jn_ /'^-r^i-r)r 



p.w-;^i:"n(-i^)f*i. (3.32) 



where is the lattice analog of the form factor (1.3) 

sin((?i - gj)/2) 



Fnh] = n 



sinh((7, + 7,)/2)^ 



Fo[q] = 1. (3.33) 



I emphasize once more that the momentum spectrum over which the summation (3.32) is 
extended (the intermediate states summing up) occurs to be boson one contrary to the 
fermion spectrum by which the initial values were defined, in particular the matrix A^x' 
(2.10). 

The function ry^ = r/(gi) that is contained in (3.32), decreases rapidly when the cylinder 
circumference grows. It follows from its definition (3.29) 

If dzhiTiz) 1 /" (ip(l — e~^(^) cosp) , , , ^, ^ 

rj{q) = : J> ^ ^ = - / ^\ . . ^Incoth AT^ p /2 . 3.34 

m J e-y(i)-z 77 J cosh7(g) - cosp ^ ^ ^ ^ 

\z\=l "-^^ 

One can find at A'"7o » 1 

^^^'> - e7(9) - 1^7Vsinh7o ' 

and r){q) —>■ ii N —>■ oo. 



4 The Paramagnetic region: K < Kc 

The correlation function is defined uniform through the determinant det A^'') in the whole 
region of parameters (including K < Kc) . But the method of its calculation is to be modified 
because the Wiener-Hopf sum equation technique, which is used, demands the Toeplitz 
matrix kernel of the appropriate form. Both the kernel and the logarithm of the kernel have 
to satisfy the Loran expansion condition. Meanwhile, the matrix kernel (2.10) is reforming 
compare to the ferromagnetic case (2.15), when K goes across Kc 

A{z) ^z-^B{z), 



where 



B{z) 



(1 



\ (1 - e-'yoz-^){l - e-T-^-i) 



(4.1) 
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so, rather the function B{z) than A{z) possesses appropriate analytical properties. The 
matrix Air) now has the following form 

= -7^ / -z-^+^'-'B{z). (4.2) 



27ri 

\z\=l 

Factorizing the kernel (4.1) 

B{z)=P{z)-Q{z-'), 

where the functions P{z) and Q{z) are analytic inside the circle \z\ < e'^° one can obtain 

nW(eT(«) - z) n^'Hi - e-^iih-^) 

Dt \ ( 70 + Ttt \ 0^9^'^ r\( -IN 0<q<7r , . 

P(.)=exp( — j n'^'(e^<"-^) ' n'»(l-e-^W.-) - ''''^ 

0<q<7r 0<9<7r 

Comparing with the factorized representation in the ferromagnetic case one can see that 
the corresponding products (3.4) and (4.3) differ from each other by one term: the factor 
(e'''" —z) is appeared in P{z) and the factor (1 — e~''''°2;~^) is disappeared in QizT^). It follows 
from (4.3) 

P(0) = Pr(0) = e-VA, g(0) = Qt(0) = 1, 

where the functions Pt{z) and Qt{z) are the same as in the ferromagnetic case (3.22), (3.23). 

The matrix (4.2) has the structure (A31) considered in the Appendix. Therefore its 
determinant can be represented according to the eqs. (A32)-(A35) by the following way 

detA^'-) =e'*('^+^)F(r), (4.4) 

oo 

nr) = E/2w(r). 

1=0 

The function h{r + 1) is given by eq. (A23). Writing it similarly to (3.6) and allowing for 
the difference between the definitions (4.3) and (3.4) for the corresponding functions P{z) 
and Q{z) one can obtain 

h{r + 1) = -(r + 1)/A + ln({{T) + ln(l - e^^""^"), (4.5) 

where the values A, ^ and are the same as in the ferromagnetic phase (3.7), (3.10) and 
(3.11). With the allowance (4.5) the eq. (4.4) can be rewritten as follows 

det Ai'^ = sinh((7o + 7,)/2)e'^WG'p(r), (4.6) 

where 

Gp{r) = 2e-^/^-(^°+^-)/2^(r) = Y.921+i{t), (4.7) 
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and for the coefficients g2i+i{r) in tfie expansion (4.7) one obtains from (A35) 

'■^'^')- N<i 

n' n {Z2i-Z2jyii 'ff (^2^-1 - ^2,-1)' h{z2iW{z2i)) 

X '^^^^^ ^ ^ . (4.8) 

n U{i-Z2iZ2j+iy U{z2i+lW{z2i'+i)) 
1=1 j=0 1=0 

The functions W{z) = P{z) / Q{z~^) and W~^{z~^) — Q{z)/P{z~^) entering the integrals 
(4.8) have the following form in this case 

W{z) = 2e-^/^-(^°+-^)/^Q^(.)sinh^(p)l^^''^— --i , (4.9) 

Ng cosh7(gj — cosp 

W-\z-') = le(^o+T.)/2g2(^) ly^W 1 _ (4_1Q) 

^ ^ 2 ^^^'N^ cosh7(g) -cosp ^ ^ 

It follows from (4.9), (4.10) that the integrand in (4.8) has only the simple poles inside the 
integration contour at the points 

^W=e-^(''«), qik) = ^, k = l,2,...,N; 

similar to the previous case. At these points, obviously, 

cosp = cosh 7(g), cosh7(p) — cosq, sinh^ 7(p) = — sin^ q. 

So, the integration contours in (4.8) may be squeezed and the coefficient g2i+i{r) is expressed 
through the sum of residua 

i^i + iJ-^V V smn7, / ^ sinh2((7, + ^j)/2) 

i<j 

where 

V2i+i[q] = Cf+i2-2'('+i)n ri (cosg2i-cosg2,)'x 

i=lj=i+l 
I l+l I 

n {cosq2i-i - cosq2j-iy[[sm^ q2i. (4.12) 

i=l j=i+l i=l 

The expression (4.12) differs from the corresponding ferromagnetic one (3.30). But it occurs 
that also in the paramagnetic case the symmetric part of the function V21+1 [q] - with respect 
to change of each variable with even subscript q2i by each variable with odd subscript q2j-i 
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(the number of the permutations is Cf'^^) - coincides with the even part of the function 
U2i+i[q] 

i<j ^ ^ 

The equahty analogous to the eq. (3.31) takes place 
where 

'-'I p 

[/— [g] = 2-(^'+^) E n-n^(^^S^)- 

W=±i] i<j \ ^ / 

Consequently the function V2i+i[q'] can be replaced by the function U2i+i[q] under the sum- 
mation in the r.h.s. of eq. (4.11). We obtain eventually 

n = 2Z + 1, 

Ur) = t^^->n{^yM (4.13) 

The functions 7^ = 7(^4) and rii = i]{qi) are defined by the eq. (3.29). 

The form of the coefficients gni'r) and of the function Fn[q] is surprisingly the same 
both the ferro- and paramagnetic cases. In the ferromagnetic case the correlation function 
expansion is extended over even n and for the paramagnetic case - over odd n 

00 

{a{r^)a{r2)) = (C • eT)e-'-/^ E 52/(r), K > K,; (4.15) 

1=0 

{a{r,)a{r,)) = sinhf^^^) (^ • eT)e-^/'^ £ ^2;+i(r), K < K,. (4.16) 

Notice that the even values of the circumference N were assumed in the factorization 
procedure for the functions P{z) and Q{z~^) (see eqs. (3.4) and (4.3) ). For the odd the 
Brillouin zone does not contain the q — tt point in the boson case. On the contrary this value 
appears in the fermion spectrum. It is not difficult to account for this detail, it is irrelevant 
for the final result. 
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5 Conclusion 



It should be observed first of all that, although the representations (4.15), (4.16) for the 
correlation function are obtained for different regions of the coupling parameter {K > Kc 
and K < Kc respectively), both of them are valid and equal to each other for all values of 
K at any finite cylinder circumference A^. For the finite the expansions (4.15), (4.16) are 
not series rather finite sums due to the form factor Fn[q] is identically zero for n > N. Only 
the thermodynamic limit N ^ oo splits the domain of their validity. 

The obtained lattice form factor representations possess one beautiful feature: they can 
be checked numerically or analytically. For small values of N the result may be computed 
in other way, for example with the help of the transfer matrix technique. In particular at 
N — 1 we have the one dimensional Ising chain. The elementary calculation gives 

(a(O)a(r)) = f. (5.1) 

One can see after simple transformations that both the representations (4.15) and (4.16) 
coincide with (5.1). 

Apart from its intrinsic interest the Ising model deserves scrutiny because of its relevance 
to nonperturbative explorations in quantum field theory. Mutual correspondence between 
field model and spin lattice system in scaling regime implies in our case the limit K — * K^, 
{N,M,r) — s> oo at {N/i, Mn,riJ,) = const. Mysterious process of thermodynamic or 
scaling limit may be observed in detail proceeding from the explicit expressions obtained in 
this work. 

The IM correlation function {a{r)a{r')) in scaling hmit corresponds to the quantum 

field model two point Green function G{x,t; L, f3). The variables of the Green function are 
assumed to be scaled on mass: m = 1, x = fi\rx — r'^\ is the spatial distance, t = fi\ry — 
is the temporal one. L = /iM is the volume, P = /iN is the temperature. The boundary 
conditions along the temporal coordinate are quite definite: periodic for the boson field, 
antiperiodic for the fermion field. This is the issue of the appropriate formulation of the 
quantum field theory at finite temperature. On the contrary, the spatial boundary condition 
may be arbitrary. In any case the presence of boundaries breaks down explicitly the Lorentz 
invariance (the rotational invariance in Euclidian metric case). The spatial and temporal 
coordinates are not equivalent in general. In this sense the configuration of correlating 
spins considered in this paper corresponds to the equal time Green function G{x,0;oo, (3) 
in the infinite volume at finite temperature. With some stipulations this function may be 
viewed also as that at zero temperature in finite volume L = (with periodicity in spatial 
coordinate). One has to keep in mind in addition that this is not equal time rather equal 
site Green function G(0, t; L, oo). It is necessary to consider the correlating spins placed on 
a cylinder circumference to compute the equal time Green function in finite volume and zero 
temperature. The corresponding problem may be solved by the method used in this work 
but with the Toeplitz determinant technique being properly modified. 

Meanwhile, the representation (4.13)-(4.16) for the IM correlation function has such 
transparent structure that the conjecture for the case of general displacement of correlating 
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spins suggests itself. The simplest one is 



{a{0)a{r^,ry)) 

It is amazing that this generalization is in exact agreement with the transfer matrix results 
for the N-rows Ising chains. We have checked it analytically for A?" = 2, 3, 4 and numerically 
for = 5, 6. 

I thank V.Shadura for the fruitful discussions and O.Lisovy for the assistance in compu- 
tations. I am grateful to Prof. B.M. McCoy for his interest to my work and useful comments. 
This work was supported by the INTAS Program (Grant INTAS-97-1312). 



g-n/A 



-rxfj-iryqj-rij 



sinh jj 



(5.2) 
(5.3) 



Appendix 

The elements of the Toeplitz matrix arranged along the diagonals are equal 

A^^j:> = Afe for k ^ X - x' ^ 0, ±1, . . . , ±r - 1 , (Al) 

\ 



Air) 





A-i 


A_2 ■ 


■ • A-r+l 




Ao 


A-i ■ 


• • A^r+2 


A2 


A, 


Ao • 


■ ■ A^r+3 



(A2) 



y Ar-i Ar-2 Ar-3 ■■■ Aq J 

The superscript in A^^) shows the matrix dimension. If there is the function of the complex 
variable A{z) analytical in the ring a < |^| < {a < 1) such that the coefficients Aj^ of 
its Loran series 

00 

Aiz) = Yl Akz'' (A3) 

k=—oo 

coincide with the matrix element (A2) then the matrix (Al) can be represented in the form 



dz 
2m J z 

\z\=l 



-x+x' 



A{z). 



(A4) 



The function A[z) in the integrand of (A4) is called the kernel of the Toeplitz matrix A^'^^ . 

If the inverse matrix is known then the calculation of the determinant reduces to calcu- 
lation of trace with help of the identity 

-1 d 



Indet A^'^HA) = JdXSp {A^'-^X))-' —A'^'-\X) 



+ const, 



(A5) 



where A is a parameter of which the matrix elements depend on explicitly. 
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The inversion procedure for the Toephtz matrix simphfics if both the kernel (A3) and 
In A{z) may be expanded in Loran series. If so, the kernel A{z) can be expressed in factorized 
form 

A{z) = Piz)-Q{z-'). (A6) 
The functions P{z) and Q{z) are analytic inside the circle \z\ < a~^: 



In P{z)^V In A{z) , In Qiz'^) ^Vln A{z) . 



(A7) 



The projection operators V and V extract the analytic part of the given function v{z) inside 
and outside the circle |-z| = 1 correspondingly 



V^=V, V^ = V, VV = VV = 0, V + V = l, 



Vv{z) 



1 f d^vjO 
2m J £ — z 



, Vv{z) 



1 f d^vjO 
2m J z — £ 

\z\=l 



(A8) 



The problem of the Toeplitz matrix inversion can be reduced to solving the system of 
the Wienner-Hopf sum equations [3] . In terms of the projection operators (A8) the matrix 
elements of the inverse matrix {A^^^)~l., are 



(A^'Y.i' = ^^i3Ty^«^-^(^)[l - Vw{z)Vw-\zr'V^^^^ 



2=0 



z=0 



Piz) 



1-V[1- w{z)Vw-\z)Vr'w{z)Vw-\z) \V-^f-^, 



z=0 



where 



w{z)^z'W{z), W{z) 



P{z) 



(A9) 
(AlO) 
, (All) 

(A12) 



Q{z-^)' 

Each of the alternative but equivalent representations (A9)-(All) for the inverse matrix 
{A^'^'>)~^^, are preferable for the calculation of different matrix elements. In particular, the 
eq. (A9) is suitable for the calculation of the element (^*^''^)r-i,0) ^Q- (^H) ~ ^r the 
element {A^''^)^}^: 



^^^^^^'^'^ = P(0)g(0) {^ - - ^i^)'Pw-\z)V]-'w{^)Vw-\z)] 



z=0 



(A13) 
(A14) 



Expanding the inverse operators that are contained in the r.h.s. of (A13) and (A14) into 
geometrical progression 



[1 - rw{z)rw-\z)r' = Y,['Pw{z)Vw-\z)]\ 

1=0 
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1=0 

and using the integral representations for the projection operators V, V (A8) one can obtain 
for (A13) with allowance (A12) 

-1 oo 

(^^'■^).--i,o = gi(^E«2m(r), (A15) 

21+1 I 

= j^-^ / 4^ -f^ . (A16) 

i=l 1=0 



\zi\<l Z1Z21+1 fl (1 - ZiZi+i) n W{z2iXi) 
i- 

In particular, the coefficients ai(r) and 03 (r) are 

dzz^-"^ 



2m J W{z-^y 

\z\=l 

, , 1 / d{ZiZ2Z^) {Z1Z2Z3Y W{Z2 

03 (r) - 



{2m)\ J ZiZ3{l - ZiZ2){l - Z2Zs) W{zi')W{z3^y 
Analogously one can obtain for the eq. (A14) 

T , nidZiZl) hw(z2i-i) 

^ ^ \zi\<i Z1Z21 n (1 - ^i^j+i) n w{z2i) 

1=1 i=l 

and, in particular, 

1 / d{ziZ2){ziZ2Y W{zi) 



1 / rf(zi2;2Z32;4)(2;i2;22;32;4)'' W{zi)W{z3) 



(27rz)2 / zM'^ - Z1Z2) W{z2y 

\zi\<l 



{2tii)\J z,z^il-ZiZ2)il-Z2Z3)il-ZsZ^) W{z2')W{z^'y 
Consider now the auxiliary Toeplitz matrix A^'^\X) depending on the parameter A, 

A,,AX)^^^ f ^z-^+^'Piz)Q\z-'), 0<A<1, (A19) 

where the functions P{z) and Q{z) are the same as in the kernel of the matrix A^'^\ It is 
seen from the definition (A19) that at A = 1 
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Zm J z 

\z\=l 

and at A = 

A.,AQ) = 7^ / ^z-^+^'P{z). (A20) 



27ri 

The matrix (A20) is the triangular one: all its elements to the right of the main diagonal 
are equal to zero 

^.,.(0) = P(0), A,,,iQ) = a x'>x. 

Therefore 

dct A(')(0) = P"(0). (A21) 

The integration constant in (A5) is defined with help of (A21). 

Let us select the contribution h{r) in In det A^'^^ which does not vanish when the matrix 
dimension grows 

In det = h{r) + Ah{r), Ah{r) ^ if r ^ oo. 

It is not difficult to observe that only the first term of the inverse operator expansion in (AlO) 
contributes in h{r). This fact results in drastic simplification of corresponding computations. 
The inverse matrix (AlO) is 



,x r 



P{z) Q\z-^) P{z) 

d{ziZ2Z3)Zi''{z2/z3)z^''' P{z2) 



{2m)\ T , Z,{Z2-Z,){Z2-ZS) P{Z^)Q\Z2^)P{Z3)' 



The derivative of the matrix (A19) is 



\z\=l 

and finally one has for the trace in (A5) 

Sp[(A(')(A))-|A<')(A)] = pi^ / Xid, 

[1 - {z^/ziYWl - {z3/z^Y]{z2/z3Y P{z2)P{z,)Q\z^^)\nQ{z^^) ^ 

{Z2 - Zi){z2 - Z3){Zi - Zi){Zi - Z^) P{ZI)Q^{Z2^)P{Z3) 

= rlnQ(0) + -^ (f dz In Q{z-^)^ In P{z) + o{r-^). (A22) 

\z\=l 



_ X 

i=l 
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It is seen from (A22) that nonvanishing contribution in the trace occurs to be independent 
on the parameter A. Therefore, with allowance (A21) we obtain 

/i(r) = r[lnP(0) + lng(0)] + -^ <f dzlnQ{z-^)-^lnP{z). (A23) 

\z\=l 

For the calculation of Ah{r) the eq. (A5) may be used by keeping the next terms in the 
expansion of (AlO). But it is more convenient for this purpose to exploit the specific property 
of the Toeplitz matrix i.e. 

^ ^0'° det A(r+i)- 

Therefore 

r+k 

detA^'-) =detA('^+*^) H {A^''^)o}>, 

s=r+l 

and after substitution (A17) we obtain 

r+fe ^ oo s 

det/lM =e'*(^+*^)+^'^(^+'=)[P(0)g(0)]-'- n (l-^a2/(s)). (A24) 

s=r+l ^ 1=1 ' 

It follows from (A23) that 

^ [p(0)Q(0)]*^e'*('"\ 

so, the factors P(0) and (5(0) in (A24) cancel. Due to this fact the limit A; — > oo in (A24) 
can be taken 

dQlA^'^ ^G{r)-Q^^'\ (A25) 

where 

n (l-E«2K«) ■ (A26) 

s=r+l ^ l=\ ' 

Expanding the product of sums in the r.h.s. of (A26) 



OC/OO s CXD OO/OC 



G{r) = \- y: e e e^^k^o e«2/(52) (A27) 

si=r+l ^1=1 ' s\=r+\ S2=si+1 ^;=1 ' \=\ ' 

summing up over Sj and collecting the terms with the definite multiplicity of integration, 
one can represent (A27) in the form 

oo 

G{r) = Y.g^,{r). (A28) 

The coefficients g2i{r) can be obtained from (A27) directly but it is more simple to compute 
them by use the following recursion relation 

I— I oo 

92iir) + J2 E 92kis)a2(i-k)is) = 0. (A29) 

fe=0 s=r+l 
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It is not difficult to deduce this relation from (A27). Using the eq. (A18) for 02/(5) and 
initial condition go{s) — 1 one can obtain from (A29) 



52/ (r) 



(-1)' 



21 i-i I I 

n {dzi o n n [(^2^-1 - ^2^-1)^(^2^ - Z2jf] n w{z2i-i) 

i=l 1=1 j=i+l i=i 



in particular, 
92{r) = 



\Zi\<l 



1 



I I 

1 n 

=1 j=i 



n n{i-z2i-iZ2,)^ 



n w{z^^] 

i=l 



, (A30) 



d{ziZ2) {Z,Z2Y W{zi] 



{2my J {l-zrZ2f W{z^') 

\zi\<l 



1 



d{ZiZ2Z^Zi) {ZiZ2Z3ZiY(zi - Z^f{z2 - Z^"^ ^^(^1)1^(^3) 



4(27rzr, / {\ - z^Z2)\\ - Z2Z^f{\ - z^z^fil - z^z^fW{z2^)W{zl^)' 



Note that the integrand of (A30) is symmetric both with respect to permutations of the even 
subscript variables and with respect to permutations odd subscript variables. 

The structure of Toeplitz matrix is altered in paramagnetic region. Here it is necessary 
to evaluate the determinant of the matrix which has the following form 







Ao 


A, ■ 


•• Ar-2 




A-2 




^0 • 








A-2 


^-1 • 


■ ■ 




^ A-r 


A-r+1 


A-r+2 • 





(A31) 



dz 

f — ' 
2Tii J z 

\z\=l 



-x+x'-l 



A{z) 



The technique stated above is not valid because the kernel of the matrix B^^^ 

B{z) = z-^A{z) 

docs not satisfy the Loran expansion condition for \nB{z) if \nA{z) docs satisfy. Never- 
theless, the problem can be reduced to that considered above. Really, it follows from the 
definition of the inverse matrix element through the determinant and cofactor 

det5(^) = (-l)'-(^('-+^))-oMet^(^+^). 

Using the eqs. (A15), (A16) for the inverse matrix and the eqs. (A25), (A28), (A30) for the 
determinant one can obtain 



detB^'^ = (-l)'-e'^(''+^)F(r), 



(A32) 
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where 

-1 oo 

F{r) = G(r + 1) • {A^^-''^);i = ^EWW, (A33) 

I 

f2i+i{r) = J2 «2fc+i(r + l)g2(i-k){r + 1). (A34) 

fe=0 

All summands in the r.h.s. of (A34) are the integrals of 2/ + 1 multiplicity. After reduction 
the integrands to common denominator and symmetrization over even subscript variables 
and over odd subscript variables one can obtain 



(-1)' f 

\Zi\<l 



u n {z2i - Z2jf n n (^2^-1 - ^2,-1)' n {z2iW{z2i)) 

^ t=l j=a+l 1=1 3=1+1 j=i /'AQc;^ 

I I I ' \^O0) 

n n (1 - ;^2^^2,+i)2 n (^2^+iW^(^2m)) 

i=l j=0 i=0 



and, in particular, 



J. , s 1 r dzz 
h[r) = — 



27ri J ziy(z-i)' 

|z|=l 



. / X _ 1 /■ 'i(^1^2^3) (^1^2^3)''(^1 - ^3)^ 2:21^(2:2) 

2{2mff {l-zrZ2ni-Z2Z^f{l-z^z^yz,z^W{z^^)W{z^^y 

\Zi\<l 
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